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FIRST-PRINCIPLES CALCULATIONS OF THE 
STRUCTURAL, ELECTRONIC AND ELASTIC PROPERTIES 
OF SrGeO3 AND SrZrO3 CUBIC PEROVSKITES 


Mikhail G. BRIK!, Nicolae M. AVRAM?, Calin N. AVRAM? 


Abstract. The structural, electronic, elastic and thermodynamic properties of two cubic 
perovskites — SrGeO3 and SrZrO3 — were calculated using the first-principles methods for 
the pressure range from 0 to 25 GPa. Comparison of the calculated results with other 
literature data (whether available) yielded good agreement. Dependencies of all calculated 
properties — such as lattice constants, relative change of the unit cell volume, elastic 
constants and Debye temperature — on pressure were obtained, which enables reliable 
estimations of all these parameters for any value of hydrostatic pressure in the studied 
range. 
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1. Introduction 


Perovskites is a common name for a large group of compounds with the general 
chemical formulas ABX3 (if A and B are the di- and tetravalent metals, then X is 
oxygen; if A and B are the mono- and divalent metals, then X is a halogen). The 


possible combinations of these chemical elements should satisfy the stability 
ratrx 


v2 (rgtrx) 
radius, the subscript i denotes an element in the chemical formula) should be 
confined in the range between 0.75 and 1 [1,2,3]. It is also possible to make the 
doubled perovskites A2BB ’'X6 by doubling a unit cell of an ABX3 perovskite along 
the c crystallographic direction [4]. Versatile chemical composition of perovskites 
and their ability to accommodate organic complexes instead of some cations may 
also lead to a variety of stable crystallographic structures, such as cubic, 
rhombohedral, tetragonal, orthorhombic, and monoclinic. 


criterion: the so called Goldschmidt tolerance factor t = (7; is the ionic 
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Applications of the perovskite compounds are numerous. Some examples that prove 
technological importance of perovskites can be mentioned as follows: photovoltaics 
[5], magnetoresistive materials [6], fuel cells [7], lighting and display materials [8], 
catalysts [9], scintillators [10,11] etc. 


Due to a very wide area of use and applications, various representatives of 
perovskite family were subjected to scrutinized theoretical studies, especially by 
means of the modern computational approaches based on the density functional 
theory (DFT). For example, surface terminations of SrTiO3 were analyzed in Ref. 
[12], systematic studies of various ABO: perovskites surfaces described in Refs. 
[13,14, 15]; bulk electronic structures for BaTiO3 and SrTiO3 were calculated in 
Refs. [16] and [17,18], respectively; bulk properties of ASnX3 (A=Rb, K; X=Cl, Br) 
were analyzed in Ref. [19] etc. Computational methods have also been proved to 
be efficient in predicting new materials with perovskite structures at extreme 
pressures [20,21]. 


Nowadays comparative studies of several or even more isostructural compounds 
acquire special importance, since they allow for uncovering certain “structure — 
property” and “property — property” trends, which can be useful guides in a search 
for new materials with improved performance. Such calculations are very popular 
now, which can be proved by a large number of publications on this topic [22,23, 
24, 25, 26, 27, 28 etc.]. In the present paper we continue this line of research and 
perform the detailed calculations of the structural, electronic, elastic and 
thermodynamic properties of two cubic perovskites SrGeO3 and SrZrO3 to see how 
the change of the second cation affects the above mentioned properties. 


In the next sections we describe the structure of the studied crystals, method of 
calculations and the obtained results before concluding a paper with a short 
summary. 


2. Crystal structure of cubic SrGeOs and SrZrOs and details of calculations 


Cubic perovskites SrGeO3 and SrZrOs crystallize in the cubic space group Pm-3m, 
No. 221 [29,30]. The Ge (or Zr) atoms are six-fold coordinated by the oxygen ions, 
which form an ideal octahedron. The Sr ions are located in the cavities between the 
SrO¢ (ZrO¢) octahedra and are 12-fold coordinated by the oxygen ions, as shown in 
Fig. la. 


All calculations were performed using the CASTEP module [31] of Materials 
Studio. Both generalized gradient approximation (GGA) [32] and local density 
approximation (LDA) [33,34] were used to treat electron correlation effects. The 
ultrasoft pseudopotentials were employed for a description of interaction between 
the ionic cores and valence electrons. The electron configurations were as follows: 
4s*4p°5s? for Sr; 4s*4p°4d°5s* for Zr; 4s*4p? for Ge and 2s?2p* for O. The plane- 
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wave bases cut-off energy was 380 eV for both crystals. The k-points grid was 
8x8x8; the convergence criteria were set as 5x10 eV/atom for energy, 0.01 eV/A 
for maximal force, 0.02 GPa for maximal stress and 5x10* A for maximal 
displacement. The calculated electronic structures were analyzed along the 
directions in the reciprocal space determined by the high-symmetry points of the 
Brillouin zone (Fig. 1b). 


To analyze the pressure effects on the calculated properties, all calculations runs 
were performed at the ambient and elevated hydrostatic pressures up to 25 GPa, 
with a step of 5 GPa. 


a) 


Fig. 1. a) One unit cell of SrZrO3. Coordination octahedra around Zr ions are shown. Drawn with 
VESTA [35]. b) Brillouin zone for a unit cell of SrZrO3 with indication of the high symmetry points, 
whose coordinates in units of the reciprocal lattice unit vectors are as follows: G(0, 0, 0); R(’4, %, 
'’); X('A, 0, 0); MC, %, 0). 


3. Results of calculations 


3.1. Structural properties 


The optimized lattice constants for the considered in the present paper perovskites 
in comparison with the experimental data and other calculations are collected in 
Table 1. 


As seen from Table 1, the SrZrOz lattice constant is greater than the SrGeO3 because 
of a large ionic radius of the six-fold coordinated Zr ion (0.72 A) in comparison 
with that for Ge (0.53 A) [36]. Agreement between the experimental and calculated 
in this work lattice constants is good; the obtained results are also consistent with 
others reported in the literature. 
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Table 1) Calculated and experimental structural constants (in A) for SrGeO3 and SrZrO3 cubic 
perovskites. 


Calc. (this work) Calc. (other works) 
GGA LDA 
SrGeO3 | 3.7946? | 3.8563 | 3.7308 | 3.8029 * 3.859 4 3.78 © 3.896! 
SrZrO3 | 4.0930° | 4.1712 | 4.0953 | 4.17942 | 4.1574 | 4.174% | 4.1575 


Exp. 


@ Ref. [29] > Ref. [30] © Ref. [37] “ Ref. [38] © Ref. [39] ‘ Ref. [40] £ Ref. [41] " Ref. [42] ‘ Ref. [43] j 
Ref. [44] 


Application of the external hydrostatic pressure will lead to a decrease of the lattice 
constants, which is illustrated by Fig. 2. 
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Fig. 2. Dependence of the calculated lattice constants on external hydrostatic pressure P for SrGeO3 
(a) and SrZrO3 (b). The calculated values are shown by symbols, the linear fits are shown by straight 
lines. The fits equations are given in the graphs. 


The lattice constants decrease is well approximated by the linear fits, with the slopes 
of about 0.005 — 0.007 A/GPa. The SrZrO3 appears to be somewhat more 
compressible than SrGeOs, which is due to the increased interionic distances and 
softer chemical bonds in the former material. 
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Fig. 3. Dependence of the calculated relative volume change of a unit cell on external hydrostatic 
pressure for SrGeO3 (a) and SrZrO3 (b). The calculated values are shown by symbols, the fits to the 
Murnaghan equation of state are shown by straight lines. The estimated bulk modulus B and its 
pressure derivative B’ are given in the graphs. 


Higher compressibility of SrZrO3 is confirmed by the data presented in Fig. 3, 
which shows the relative change of the unit cell volume V/Vo (where V and Vo are 
the unit cell volume at pressure P and ambient pressure, respectively). The values 
of the V/Vo ratio were fitted to the Murnaghan equation of state [45] 

1 


La(ie0e)™ e 


0 


where B is the bulk modulus and B’ is its pressure derivative. The obtained values 
of these elastic constants are given in Figs 3a and 3b. In addition, the experimental 
data on change of the relative volume of a unit cell for SrGeOs at different pressures 
[40] are shown in Fig. 3a by the stars. Very good agreement between the calculated 
and experimental data was demonstrated. Unfortunately, no experimental data on 
compressibility of SrZrO3 were found in the literature. 


3.2. Electronic properties 


The calculated band structures of both perovskites are shown in Fig. 4. It is known 
that the DFT-calculated band gaps are underestimated if compared with the 
experimental ones. 
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Fig. 4. Calculated electronic band structure for SrGeO3 (a) and SrZrO3 (b). 


The calculated band gap values for SrZrO3 (3.328 eV for the GGA calculations, and 
3.313 eV for the LDA calculations) are consistent with others reported in the 
literature: 3.413 eV [42], 3.31 [43]. The calculated band values for SrGeO3 (~0 eV 
for the GGA calculations, and 0.789 eV for the LDA calculations) are smaller than 
the reported value of 2.84 eV in Ref. [39], but it should be emphasized that the 
authors of Ref. [39] used the hybrid DFT method, which improves agreement with 
the experimental band gap values. The lowest states in the conduction band are 
highly dispersive, whereas the states from the valence band top are nearly flat, at 
least in the vicinity of the R and G points. 


The band gaps are of indirect character for both studied materials. The origin and 
composition of the calculated electronic bands becomes clear after a careful 
examination of the density of states (DOS) diagrams from Fig. 5. 
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Fig. 5. Calculated density of states diagrams for SrGeO3 (a) and SrZrO3 (b). 


The valence bands in both compounds are mainly composed of the oxygen 2p states. 
The conduction band in SrGeO3 is made by the Sr 5s states with a minor 
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contribution from the Ge 4s, 4p states. The conduction band in SrZrO3 consists of 
the Sr 5s states and Zr 4d, 5s states. The Sr 4s, 4p states produce sharply peaked 
maxima at about -33 eV and -15 eV, respectively. The O 2s states are located at 
about -17 eV in SrGeOs3 and at about -15 eV in SrZrO3. The Zr 4s and 4p states are 
located deeply at -47 and -25 eV in SrZrO3, whereas the Ge 4s, 4p states are at much 
higher energies in SrGeO3. 


The chemical bonding properties can be analyzed by considering the effective 
Mulliken charges and Mulliken bond orders. Because of formation of chemical 
bonds between ions in solids, the effective charges differ from those ones prescribed 
by their formal valences. Thus, the calculated charges of all ions in SrGeOs3 are (in 
units of the proton charge, GGA/LDA results are listed): Sr +1.01/+0.97; Ge 
+1.72/+1.70; O -0.91/-0.89. The same values for SrZrO3 are: Sr +1.40/+1.36; Zr 
+1.07/+1.01; O -0.82/-0.79. The charges of the Sr ions are close to the formal value 
of +2, which indicates a mostly ionic Sr — O bonding, whereas the Ge — O and Zr — 
O bonds are more covalent in comparison with the Sr — O case. The larger value of 
the bond order corresponds to the covalent character of the considered bond, 
whereas its lower values indicate the ionic nature of the bond [*°]. Thus, the Zr — O 
bond orders in SrZrO3 are 0.86/0.84, whereas the Ge — O bond orders in SrGeO3 
are 0.86/0.88, indicating approximately comparable degree of covalency of those 
bonds. The values of the Sr — O bonds are 0.13/0.14 in SrZrO3 and 0.18/0.18 in 
SrGeOs, which is another indicator of highly ionic bonding in these cases. 


3.3. Elastic properties 


In addition to the extraction of the bulk moduli from the Murnaghan equation of 
state, direct calculations of the elastic constants were performed at the ambient and 
elevated pressures. The elastic properties of a solid are described by a 6 x 6 matrix 
of its elastic constants Cj. In the case of a cubic crystalline solid, only three elastic 
constants are needed: C11, Ci2, C44, with the following relations: C1 = C22 = C33, 
Ci2 = C23 = C31, C44 = C55 = Coo, and all remaining constants not listed here are 
equal to zero [47]. 


The calculated elastic constants values at the ambient pressure are collected in Table 
Zs 
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Table 2) Calculated and experimental elastic constants (in GPa) for SrGeO3 and SrZrO3 cubic 
perovskites. 


Calc. (this work, GGA/LDA) Calc. (other works) 
Cu C12 Cag B Cu Cx Cag B 
SrGeO3 | 284/373 | 96/130 | 138/130 | 158/211 - - - - 
310°, da 74°, 151°, 
309°, 71°, 74, 155°, 
SrZrO3 | 309/367 | 73/76 74/76 | 152/173 307°, 68°, TAS, 148°, 
3399 719 774 1604 


@ Ref. [41] ® Ref. [42] © Ref. [44] ¢ Ref. [48] 


No other data for SrGeO3 were found in the literature, but high consistency between 
our calculated values and other results for SrZrO3 allows to consider the SrGeO3 
values to be reliable. 


It is worthwhile mentioning that the values of the bulk moduli B derived from the 
“volume — pressure” dependence in Fig. 3 are in good agreement with the data from 
Table 2. 


Fig. 6 shows the behavior of all calculated lattice constants with increased pressure. 
The calculated values are shown by the symbols (solid ones for the GGA values, 
open ones for the LDA values). All data were linearly fitted; the fit lines are shown 
as solid for the GGA data and dashed for the LDA data. The linear fits equations 
are collected in Table 3. 
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Fig. 6. Calculated elastic constants as the functions of external hydrostatic pressure for SrGeO3 (a) 
and SrZrO3 (b). See the text for more details. 
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Table 3) Dependences of elastic constants on pressure P (in GPa) for SrGeO3 and SrZrO3 cubic 


perovskites (equations of straight lines from Fig. 6). 


SrGeO3 SrZrO3 
GGA LDA GGA LDA 
Cu 286.9 + 8.0 P 374.9+7.8P 311.0+8.1P 368.3 + 8.4 P 
Cr 97.5 + 3.0 P 130.9 + 2.8 P 73.0+1.9P 76.7+2.1 P 
C4 139.44+ 1.6P 175.5+1.5P 74.4 +0.3 P 76.2 +0.3 P 
B 160.2 + 4.6 P 212.24+4.5 P 152.5+4.0 P 173.8 +4.2 P 


Equations given in Table 3 allow for a reliable estimation of the elastic constants at 
any value of the external hydrostatic pressure in the range from 0 to 25 GPa. 


3.4. Thermodynamic properties 


After the values of the elastic constants for a given solid are determined, it is 
possible to estimate the Debye temperature Op, which is one of very important 
parameters of solids, proportional to the highest energy of the acoustic phonons. 
The following equations can be used [47,49,50]: 


aah a, i 


kl4n\ M 


where / and k are the Planck’s and Boltzmann’s constants, Na is the Avogadro’s 
number, p stands for the crystal’s density, M equals to the molecular weight, and n 
is the number of atoms per one formula unit (five in this case). The last entry, vm, is 
the mean sound velocity, which can be calculated if the longitudinal v; and 
transverse v; sound velocities are known: 


Ric (3) 
1(2 1 
Ym = Aaa) 
3B+4G G 4 
vy = 3p , Ut = ff 2 ( 
B is the bulk modulus, G = ne is the isotropic shear modulus, which is an average of the 


Voight’s shear modulus Gy and the Reuss’s shear modulus Gr. These parameters are 
expressed in terms of the elastic constants Cjj as: 


: 5 ’GRCua-C az Cy” 


Gy (5) 
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Application of Eqs. (2) — (5) to the calculated data on the elastic constants resulted 
in the values of the densities, sound velocities and Debye temperatures collected in 
Tables 4-5. 


Table 4) Calculated density p (kg/m3), sound velocities v;, v;, vm (m/s) Debye temperatures @p (in 
K) for different hydrostatic pressures P (in GPa) in SrGeO3 cubic perovskite. 


GGA LDA 
p Vr v1 Vn Op p Vt vi Vin Op 
O | 6028 | 4430 | 7237 | 4891 646 6657 | 4766 | 7873 | 5268 719 
5 6206 | 4585 | 7605 | 5069 676 6807 | 4878 | 8140 | 5397 742 
10 | 6365 | 4697 | 7870 | 5199 699 6946 | 4976 | 8367 | 5509 763 
15 | 6511 | 4809 | 8143 | 5328 722 7077 _ | 5090 | 8609 | 5639 786 
20 | 6647 | 4904 | 8383 | 5438 742 7199 | 5169 | 8808 | 5730 803 
25 | 6774 | 5002 | 8591 5549 762 7315 | 5243 | 8982 | 5815 819 


P 


Table 5) Calculated density p (kg/m?), sound velocities v,, v7, Vm (m/s) Debye temperatures Op (in 
K) for different hydrostatic pressures P (in GPa) in SrZrO3 cubic perovskite. 


GGA LDA 

p vr vi Vn Op p Vv; VI Vin Op 
0 5188 4148 7228 4608 563 5482 4247 7457 4720 587 
5 | 5351 | 4232 | 7485 | 4707 581 5631 | 4315 | 7734 | 4805 603 
10 | 5499 | 4303 | 7720 | 4791 596 5769 | 4369 | 7908 | 4868 616 
15 | 5636 | 4355 | 7954 | 4856 609 5898 | 4415 | 8145 | 4927 628 
20 | 5763 | 4405 | 8153 | 4916 622 6018 | 4478 | 8323 | 5000 641 
25 | 5883 | 4461 | 8308 | 4982 634 6132 | 4520 | 8538 | 5052 652 


P 


As seen from Tables 4-5, all calculated parameters increase with increased 
pressures. The Debye temperature increases linearly with pressure, as is evidenced 
by Fig. 7. 
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Fig. 7. Calculated Debye temperatures as the functions of external hydrostatic pressure for SrGeO3 
(a) and SrZrO3 (b). The calculated values are shown by symbols, the linear fits are shown by straight 
lines. The fits equations are given in the graphs. 
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4. Conclusions 


Detailed first-principles analysis of the pressure effects on the structural, electronic, 
elastic and thermodynamic properties of two cubic perovskites — SrGeO3 and 
SrZrO3 — was performed in the present paper. All calculated values — lattice 
constants, relative change of a unit cell volume, elastic constants, and Debye 
temperature — were fitted to the linear functions of pressure. The obtained 
dependencies allow to estimate all these quantities for any value of pressure in the 
studied range from 0 to 25 GPa. 
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